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Abstract 

We prove the existence of a large class of dynamical solutions to the 
Einstein-Euler equations for which the fluid density and spatial three- 
velocity converge to a solution of the Poisson-Euler equations of Newto- 
nian gravity. The results presented here generalize those of [10] to allow 
for a larger class of initial data. As in [10], the proof is based on a non- 
local symmetric hyperbolic formulation of the Einstein-Euler equations 
which contain a singular parameter e = vt / c with vt a characteristic 
speed associated to the fluid and c the speed of light. Energy and disper- 
sive estimates on weighted Sobolev spaces are the main technical tools 
used to analyze the solutions in the singular limit e \j 0. 



1 Introduction 



The Einstein-Euler equations, which govern a gravitating perfect fluid, are 
given by 

where 

= {p + c~^pyv^ +pg'^, 

with p the fluid density, p the fluid pressure, f * the fluid four- velocity normal- 
ized by v'^Vi = —cP, c the speed of light, and G the Newtonian gravitational 
constant. Defining 

Vt 

e = — , 
c 

where vt is a typical speed associated with the fluid, the Einstein-Euler 
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equations, upon suitable rescaling [10], can be written in the form 

(1) G'^ = 2t^T^ and ViT^ = 0, 
where 

T*-' = {p + e^pWv^ +pg^^ and v^Vi = 

In this formulation, the fluid four- velocity v^, the fluid density p, the fluid 
pressure p, the metric gij, and the coordinates (x*) i = 1,...,4 are di- 
mensionless. By assumption, the (x*) are global Cartesian coordinates on 
spacetime M = M'^ x [0, T), where the (x^) {I = 1,2,3) are spatial coor- 
dinates that cover M^, and t = x'^/vt is a Newtonian time coordinate that 
covers the interval [0,T). By a choice of units, we set vt = 

The Newtonian limit for the Einstein-Euler system refers to the limit of 
solutions of the Einstein-Euler systems in the limit e \ 0. In this limit, 
one expects that under reasonable assumptions solutions of the Einstein- 
Euler system should converge to a solution of the Poisson-Euler equations 
of Newtonian gravity: 

(2) dtp + di{pw') = {di:=d,i), 

(3) pidtw'' + w^diw^) = -{pd-^^ + d-^p) {d'^ := S'^^di), 

(4) A$ = /5 (A:=5"a7aj), 

where p, p, and w"^ are the fluid density, pressure, and three-velocity, re- 
spectively. 

The difficulty of analyzing the Newtonian limit arises from the fact that 
the limit e \ is singular. The first general rigorous result on the Newtonian 
limit without any symmetry assumptions is [13]. There, it is shown that 
there exists a wide class of solutions to the Einstein equations coupled to 
Vlasov matter that have a well defined Newtonian limit as e \ 0. This 
work is based on an elliptic-hyperbolic formulation of the Einstein- Vlasov 
equations in a maximal slicing gauge. In [10], we used a different approach 
to prove existence of a large class of non-stationary solutions to the Einstein- 
Euler equations which have a Newtonian limit. 

The main aim of this article is to establish the existence of a Newtonian 
limit for solutions to the Einstein-Euler equations under weaker conditions 
on the initial data as compared to either [10] or [13]. The purpose for this is 
twofold. First, it is of both theoretical and practical interest to understand 
the most general situations possible for which Newtonian gravity provides 
an acceptable approximation to full Einstein theory. Second, the techniques 
developed here can be used to improve the results of [11] on the existence of 
post-Newtonian expansions. In [11], it was shown that there exists a class of 
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solutions to the Einstein-Euler equations that have a first post-Newtonian 
expansion. Using the methods here, this can be improved to the second 
post-Newtonian order. We will report on this in a separate article. 

In this article, we follow the approach of [10] to analyze the limit e \j 
of solutions to the Einstein-Euler equations. This requires that we replace 
the metric gij and fluid velocity with new variables that are compatible 
with the limit e \ 0. The new gravitational variable is a density u*-' defined 
via the formula 

(5) a'' = , ' 



where 



V-detCQ) 



(0) „ „ „ -1 +'''nu"« 







From these formulas, it not difficult to see that the density u*-' is equivalent 
to the metric gij for e > 0, and is well defined at e = 0. For the fluid, a new 
velocity variable w'^ is deflned by 

(7) = and w'^ = . 

For technical reasons, we assume an isentropic equation of state 

(8) p = 

for the fluid where K G M>o, n G N. This allows us to use a technique of 
Makino [9] to regularize the fluid equations by the use of the fluid density 
variable a deflned by 

^ " (4Kn(n + 

The resulting system can be put into a symmetric hyperbolic system that 
is regular across the fluid-vacuum interface. In this way, it is possible to 
construct solutions to the Einstein-Euler equations that represent compact 
gravitating fluid bodies (i.e. stars) both in the Newtonian and relativistic 
setting [9, 12]. 

The main point of introducing the gravitational- matter variables {u*-' ,w^,a} 
is that in a harmonic gauge the Einstein-Euler equations can be cast into a 
singular (non-local) symmetric hyperbolic system of the form 

(10) b°{eW)dtW = ^c^diW + b^{e, W)diW + F(e, W). 

As shown in [10], for appropriately chosen initial data this form is suitable 
to derive e-independent energy estimates that can be use to analyze the 
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behavior of the solutions as e \ 0, and extract a Newtonian Umit. We 
note that singular hyperbolic systems of the form (10) have been extensively 
studied [3, 6, 7, 14, 15], but, as discussed in [10], are not directly applicable 
to the Einstein-Euler equations due to initial data that does not lie in the 
standard Sobolev space H^{M.^). 

For general initial data, the e-independent energy estimates from [10] 
are not enough to control the solution in the limit e \j 0. In this paper, we 
show that when the energy estimates are used in conjunction with dispersive 
estimates for the wave equation a larger class of initial data can be chosen 
so that the resulting solutions still have a Newtonian limit. However, unlike 
the situation in [10], the gravitational variables do not converge to an e- 
independent limit. Instead, they converge to a solution of a singular e- 
dependent wave equation. Following the terminology used in other singular 
hyperbolic problems [15], we refer to this type of limit as a fast limit. We 
note that dispersive wave estimates have been used previously in a similar 
fashion to analyze the (singular) incompressible limit for the Euler equations 
[5, 16]. 

The precise statement of the existence of a Newtonian limit is contained in 
the following Theorem which is the main result of this article. A proof can 
be found in section 6. A definition of the weighted spaces H^^ [H^ := 
H^i) can be found in Appendix A of [10]. We also define Xts ^^^ := 

n£^c^([o,T),F^^). 

Theorem 1.1. Suppose -I < 5 < -1/2, s £ Z>2, R > 0, k € I^ys+s, 
a,w^ G He,, supp a C Br, i^'^ G H'^^^ , Ia^ ^ H'^^,. Then there exists a 

o o o 

T > 0, eo > 0, and maps 

unt) : u'Jit) - u'Ho), div^iit), mnt) e XT,s,k,5^i O < e < eo, 
^l^it) : ~x^^{t) - u^^(0), dfx^^it), dtK^{t) G Xt,sAS-i < e < eo, 
p,{t), w\{t) e XT^s,k,5~i < e < eo, 

p{t), w\t) e XT.s,k,5^1, ^{t) G XT,s,k+2,S with dt^{t) E XT,s,k+l,5-l, 
such that 

(i) the triple {u^ {x,t), {x,t), w\ {x,t)] determines a solution to the Einstein- 
Euler equations (1) in the harmonic gauge for < e < eo on the 
spacetime region {x^ ,t = x^) € D = x [0, T) with ADM mass given 
by 

"T-ADM = / p{x,0)dx^ + 0{e^), 
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(ii) {p{x,t),w^ {x,t),^{x,t)} is a solution to the Poisson-Euler equations 
(2)-(4) with initial data p\t=o = {4:Kn{n + l))~^a^'", w^\t=o = w\ 

o o 

(iii) Ue{x,t) is a solution to the wave equation 

e^d^iii^ - AUl^ = -5l6{~p + eHlSld^^, 
with initial conditions 

(iv) 

IIPe(t) - p{t)\\m~2 + \\wi{t) - w\t)\\H,-2 + \\wt{t)\\H^-2 < e, 
and 

WK'it) - S:^IIl6 + Wdiunt) - dj^l^t)\\j,,^, + Wedtunt) - edt^nt)\\H.-2 < e 
for all {t,e) G [0,T) x (0,eo]. 

From the above theorem, the interpretation of the hmiting solution is 
clear. The {p,w^} satisfies the standard Poisson-Euler equations of Newto- 
nian gravity with the obvious interpretation as the fluid density and three- 
velocity, while the represent high frequency gravitational radiation prop- 
agating on a flat background with the fluid density and Newtonian potential 
acting as source terms. 



2 Reduced Einstein Equations 



To aid in deriving the appropriate symmetric hyperbolic system for the 
gravitational variables, we temporarily introduce a new set of coordinates 
related to old ones by the simple rescaling 

X — X ^ X — X j 

and let 

d^ = — d, = — 

dx^ ' dx^ 

In the new coordinates, the metric cjij and its inverse g'^^ are given by 

Next, we consider the metric density 

(12) Q'^ = ^A^\9'^ where \g\ = - det{gij) . 
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We note that the metric g^^ is related to the density g*-' by the following 
formula 

(13) g'^ = ^=t^ where |5| = -det0'^ 



and hence 



(14) {g'- 



1 rg'J 60^4 



\g\ V^B ^ 



To obtain a gravitational variable that is regular and non-trivial in the limit 
e \ 0, we define 

(15) u*J:=_L(g*._^*.) 
where 

'I3X3 



VI ) I -1^ 

is the Minkowski metric density. As stated in the introduction, for e > 0, 
the metric gij can be recovered from the density u*-' via the formulas (5)-(6). 
In the (x*) coordinate system, the Christofell symbols are given by 

(16) flj = [t'"{2-QifQ,p - -QifQ^p)dmu''' + mpS^^dj^u'-P - 20,(ia,)U^^)) . 

These are related to the Christofell symbols in the (x*) coordinate system 
as follows 

-pA _ ^-2f,A -p4 _ -lf^4 -p4 _ f,4 
(,i g i 44 — e i 44 , i 44 — e i 44 , i ^4 — i ^4 , 

(18) T% = eT\B. r^4 = e-^f^4 and T^c = 



Using (16), a straightforward calculation shows that the Einstein tensor 
G*-^ is given in terms of the density u^^ by 

(19) := ^\g\ G'^ = Q^^BlfU'^ + [A^^ + B^^' + + I)^^' 
where 

(20) |5| = -det(r^), 

(21) A'^ = 2(l0fc,0^„ - -QkmQln) {n'" - irr^)5pU*^'5,U™" , 

(22) B'^ = 40fc,(2r('5™u^'^^4u''" - 5^u^"9nu™' - 0™"a^ii^'=4u^O ' 

(23) C'^ = 'i{dku'^deu''^ - Bku'^Beu^'') , 

(24) Z5*-' :=0'^a|,u'^^-24V'^V^'- 
To fix the gauge, we assume that 

(25) Biu'^ = 0. 
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For e > 0, this is equivalent to to the harmonic gauge 

(26) diQ'^ = di{y^-detigM)g'^) = 0. 

Setting 

(27) := g'^ - D'^ = Q'^^dlfU'i + [A'^ + B'i + 
and 

the Einstein equations G*-' = 2e^T*-' in the gauge (25) become 

(28) Qt='T''. 

We win refer to these as the reduced Einstein equations. 

To write the reduced Einstein equations in first order form, we introduce 
the variables 

' diu'^ ifk = I 
ediU^^ ifk = A ■ 
The reduced Einstein equations then become 

-f^xC^ = 2t'dixC^ + -g'-^diiCi + e^{A'^ + B'^ + - V' , 
or equivalently 



u/ := Buu'^ 



1. 

"a 



e 

Next, we define 

(29) u^^:=e#, u^= u^, 
and let 

V = { {r'^) G M4X4I det(r/'^' + 4r^-') > } . 
Then using vector notation 

the reduced Einstein equations take the form 

(30) A\eu)d4u'^ = -C^dju'^ + A^{u)diu'^ + F'^ie, u) - -(r^ 0, 0)^, 
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where 

'1 - 4eu44 0> 

(31) A^{eu) = \ (^^-^ + 4eu^'^ | , 

1, 







^ 6^-^ 0\ 


(32) 


C'= j 


, 






yO 0/ 






/8u^^ 4u^'^ O' 


(33) 


A\u) = 1 


4u^-^ 






V 


(34) 







and 

(35) F^^' (u, eu) = (A^-'' + B'^ + &^ , 0, 0)^. 

Here we are using the notation 

u = (u*-' ) and Ufc = (%^) . 

The stress-energy tensor is given in terms of the u variable by 

,..,.,,v).-i^(7 2).-i.(-:^2) 

which we can write as 

(37) 7(^')=(2 .-v)+^' 

where 

^'^ )-P[\-g\v'v^ e-^[i\-g\-l){vr + {ivr-i)]J^ 
(38) 

|-| A/O + e^vWv^ + |g|"^/^p((5^"^ + 4eu^'^) epv^v"^ + 4e|5|"^/2pu^^ 
^1^1 1^ + 4e|5|-V2pu^-^ p(i;^)2 + \g\-^''^p{-l + 4eu44^ 

Letting (see (7)) 
(39) w = (q,w;*)^, 
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we can decompose 5*-' as 

(40) S'^ = S'o^ + eSi\ 

where 

5o^(u, w,eu,ew) = 

(A^\ ( \g\w\l + ew^) \ 

^ ' ^\\g\w\l + ew^) e-\{\g\-l){l + ew^f + {{l + ew^f -I)]) ' 

and 

5{-'(w, eu, ew) = 
(42) 

l_l fpw^w^ + pew^ew^ + \g\~^^'^pQ^'^ ^£-^^(1 + ew^) + 4|5|~^/^peu^^ 
V Pew-^i^ + ew^) + 4|5|-VVu^^ p{l + ew^f + \g\-^/^pi-l + 4:eu^^] 

3 Regularized Euler equations 

In the coordinates (x*), the Euler equations are given by 

(43) ViT^ = 

where T*-' = [p + e^p)v'^v^ + pg'^^ and the fluid velocity is normalized 
according to 

(44) v,v' = -\. 

To write (43) as a symmetric hyperbolic system, we follow [2] and differen- 
tiate (44) to get 

(45) ViVjv' = and v^v^/jv' = . 
Writing out (43) explicitly, we have 

(46) {Bip + e^dip)v'v^ + {p + e^p){v^Viv' + v'ViV^) + g'^dip = . 
The operator 

projects into the subspace orthogonal to the fluid velocity v^, i.e. LjLl, = Ll, 
and L-'-v^ = 0. Using to project the Euler equations (46) into components 
parallel and orthogonal to yields, after using the relations (44)-(45), the 
following system 

(47) v'd,p+{p + e^p)L]V,v^ = 0, 

(48) MijV^VkV^ + —^L%p = , 

p + e^p ■' 
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where 



Mij = gij + 2e^ViVj . 



As discussed in the introduction, we use a Makino density variable a 
(see (9) ) to regularize the fluid equations in regions where the density and 
pressure vanish. After multiplying (47) by the square of the function 

a short calculation shows that the Makino density a and the fluid four- 
velocity v"^ satisfy 

(49) h'^v'dia + gLjViU^' = , 

(50) MijV^VkV^ + qL{dja = , 
where 

2 dp 1 2 

dp 4n^ 

is the square of the speed of sound, and 

1 

q = — — — r-a . 

2nh[ea) 

Instead of solving (49)- (50), we consider the following modified system 

(51) h^v'Bia + qLiViV^ = 0, 

(52) Mijv'^VkV^ + qL{dja + (xar " l)M^,f{^v''v^ = 0, 
Here we are using 

Xx{x) := x{x/\) A > 0, 
where x ^ C°°(M^) is a smooth cutoff function satisfying x(^) = 1 for 
< 1, x{x) = for |x| > 2, and < x{x) < 1 for ah x € M^. 

Since = and w'^ = v'^ — 1/e, we can write (51) and (52) as 

(53) a^(94W = a^djw + b 
where 

(54) 



h\l + €W^) eqL] 



and 
(56) 
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From (13), (15), (29), and (40), we find that 
(57) gij = Vij + fiji^u) 

where the fij{y) are analytic and satisfy /ij(y) = 0(|y|) as y — )• 0, while 

(16) shows that 

(58) 

^% = e (2r/i^r/j> - %%p)eu^ + 2(%p5f.eu^5 - 2%(ieujf)] + e4(eu,eu^) 

for functions /j^ (eu, eUm) that are analytic for eu € V, linear in the eUm, and 
satisfy /j^ (0, y) = 0. The expansion (57) allows us to write 



(59) 
and 
(60) 



l{ = si + e'^QikV^v^ 



for functions £:?(eu, e«;^) and 'mij{eu, ew'^) that satisfy £j{0,0) = ?Tijj(0,0) = 
0, and are analytic for eu G V. Using (57)-(60), we can express the a* and b 
as 



(61) 
(62) 

and 



-^Sj -5ijW 



'j^ + w^a{eu, ew) + aa\eu, ew), 



(63) 







ip 



2m^i) 



+ 



/ a6i(eu,ew) • eu^ 

VX4R^2( 



eu,ew) • Ufc 



We observe that the matrices a^, a, and are symmetric, and the maps 
a^, a, d^, 6i, and 62 are analytic (for eu G V) and satisfy a^(0,0) = 0, 
a^(0,0) = 0, a(0,0) = 0, 61 (0,0) = 0, and 62(0,0) = 0. This shows that 
the system (53) is symmetric hyperbolic on a region where (eu, ew) is small 
enough to ensure that is positive definite. This can always be arranged 
by taking e small enough and since we are interested in the limit e \ no 
generality is lost by assuming this. 
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4 Uniform local existence 



The combined systems (36) and (53) can be written as 
b°{eV, e^U)dtV =^c^diV + h\V, eU, eV, e^U)diV+ 

(64) fo{V, eU, eV, e^U) + eh{V, eU, eV, e^U) + -^g{V), 
where 

(65) ?7 = (0,0,u^^0,0)'^, u'^=u'^\ 

o o ' '-"'^ 

(66) V = {u^iyj,du'\a, wY , Su'^ = u'^ - eu'\ 

(08) I 

(69) 6'(V,£C/,£l/,£2(/)=f'''W „ " , 

^ ' \ ! J ! ^ Y (w, eu, ew^ 

(70) h{V, eU, eV, e^U) = f ^^L" ?' '"^^"l . 
^ ^ ^ ^ V &(u,w,eu, ew) y 

(71) /i(y, 6C/, 6y, e^c/) = - ^^)^ 

and 

(72) ff(y) = (-<^l<5ip(a),0,...,0)^. 
For initial data, we will often use the notation 

z = z\t=o ■ 



In addition to solving these evolution equations, we must also solve the fol- 
lowing constraint equations on the initial hypersurface S = {(a;^,0) | (x^) S 
M^} to get a full solution to the Einstein-Euler equations: 

(73) C-' := G^^ — T^* = (gravitational constraint equations), 

(74) W := diU^-^ = (harmonic gauge condition), 

and 

(75) Af := eviV^ H — = (fluid velocity normalization). 

e 
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To fix a region on which the system where both the evolution (64) and 
constraint equations (73)-(75) are weh defined, we note from (30), (61), 
and the invertibility of the Lorentz metric (??*-') that there exists a constant 
Kq > such that 

(76) - det(V^' + 4eu*^) > 1/16 , I + ew^ > 1/16 , 

(77) A\eu)>^l , a\eu,ew)>^l, 

Id 16 

and 

(78) \A^i(^u)\ < 16, |a^(eu,ew)| < 16 

for all |eu| < 2Kq, \ew'''\ < 2Ko, \ea\ < 2Kq. The choice of the bounds 1/16 
and 16 is somewhat arbitrary, and they can be replaced by any number of 
the form 1 /M and M for any M > 1 without changing any of the arguments 
presented in the following sections. However, since we are interested in the 
limit e \j 0, we lose nothing by assuming M = 16. 



4.1 Newtonian initial data 



To generate a one parameter family of solutions to the constraint equations 
(73)-(75) that is regular in the limit e \ 0, we use a slight variation of 
the method used in [10], which is based on previous work by Lottermoser 
[8]. Before we state the theorem, we note from (9), (8), and the weighted 
multiplication inequality (see [10] Lemma A. 8 ) that if q G (5 < 0, A; > 
3/2) then p,pe H^. 

Proposition 4.1. Suppose —l<6<0,k>3/2+l,R>0 and {p^p^uj^ ,^1"^ 
{H^Z2Y X X H^Zl X Br{H^) . Then there exists an eo > 0, an open 

neighborhood U of ip,p,w^ yji"^ and analytic maps (— eo,eo) x [/ — )• 

: {e,p,p,w^,ii-^,i^-^)^w'^,{-eo,eo)xU -^HI : {e, p,p,w^ 
4>, (—60, eo) X [/ — ^ Hg : {e, p,p,w^ ,^1'^ ,}^'^) i-> tv^ such that for each 

{p,p,w^ ,1^1 il^'') £ U, {e, p,p,w^ ,\xl ,d4U^^) is a solution to the three 
constraints 



C^' = , W = 0, and M = 0, 



where 



XX)^ d> 
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and 



e €544 

Moreover, if we let (pQ = (/>|e=o, tDg = tt)^|e=o, and Wq = w^\^=o, then 4>o, XOq, 
and Wq satisfy the equations 

A</.o = p + , An>i = -dL5i-\ and = , 

respectively. 

Proof. The proof follows from a simple adaptation of the proof of Proposition 
5.1 in [10]. □ 

Corollary 4.2. For —1 < S < —1/2, the ADM mass of the 1-parameter 
family of initial data constructed in Proposition 4.1 satisfies 



iriADM 



[ pdx^ + 0{e). 
Jr3 



Proof For fixed (p, p, , i^-^) G {H^^Zl? x hI_^ x H^zl x H^, it follows 

from Proposition 4.1 that for eo small enough, the maps 

(79) [0,eo) Be^g^jGH^ and [0,€o) B e ^ edtm, e H^Z^ 

are analytic (see (11) and (15)). Moreover, a short calculation shows that 

(80) glj = r]ij + 2e2 {r]kiu''%j - 2r]iku''^r]ij) + O(e^), 

and 

(81) edtgij = 2e^ {rjuufrnj - 2ij,kufvej) + Oie^), 



where 

(82) (u^^) 

(83) (v^i) = ( 



3" 



and 

(84) Act>o = p + dhl'-^ , Att)^ = -0^4'' ■ 

Since —1 < d < —1/2, it follows from Proposition 4.5 of [1] that the total 
ADM energy-momentum P = (P-,) for the initial data {gij, edtgij) on the 
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initial hypersurface S = {(a:;^,0)| (x^) G M^} can be calculated using the 
standard formulas 

(85) ^^ = -\j {d'gij-5'''djgiK)dS\ 

(86) ^' = \'f {KK^iJ-Kij)dS\ 
where the extrinsic curvature Kjj is given by (see (16)) 

(87) Kij = -\l^^h- 
Furthermore, the map 

(88) [0,eo) 3 (IPi) GM^ 
is smooth by (79) and Theorem 5.1 of [1]. 

By (16), (80), (81), and (87), we find the following expansions for the 
extrinsic curvature: 

(89) Kij = e2 (2uf + 6ij{uf - 6^4') + 2dju^' + 2diu^j) + 0(6^). 

The smoothness of the map (88) and the two expansions (80) and (89) show 
that the ADM energy-momentum can be expanded as 

(90) P4 = -e'/ {d'(t>o-dn'-^)dSi + 0{e^), 

(91) ¥i = e^i {SjdKtv^ -d-^wi-ditot'^ -d-^)dSj + 0{e^). 

Using the divergence theorem and (84), the energy-momentum expansions 
(90) and (91) simplify to 

(92) P4 = -e^ / pdx^ + 0(e3) and P/ = 0{e^). 

The proof of the corollary now follows from the above expansions and the 
definition 

1 



of the ADM mass. □ 
4.2 Uniform existence 



To prove local existence of solutions to (64) on a uniform time interval 
independent of e, we take the same approach as in [11] and use a non-local 
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modification of (64). The modified system is constructed as follows. First, 
we replace g(V) in (64) with 

(93) g{V) = {-6i6ixRp{a),0,...,0f , 
and we define the Newtonian potential by 

(94) A$ = XrP. 

Next, we use the Newtonian potential to define a new combined gravitational- 
matter variable W via the formula 

(95) W = V -d^, 
where 

(96) d$ = (0,(5i(5iaj$(a), 0,0,0). 

Notice that the transformation (95) leaves the matter variables unaffected. 
Consequently, we can define W by 



..i\T 



W 



and treat $ and as a function of VF. To formulate the evolution equation 
entirely in terms of VF, we need the time derivative of the $ map. So we 
define 

^W, eU, eW, e^U) = A'' (^^^^!!^n(a^(eu, ew)"! [a^(w, eu, ew)a/w 



(97) +6(u,w,eu,ew)]) 

where Il{{a,w^)'^) = a is a constant projection map. By construction, 
$ = dt^ when evaluated on a solution of (64). To fit with the above 
notation, we also define 

d^ = {0, 6i5idi^, 0, 0,0f . 

Noting that 
(98) 

b°{eV, e^U) = 6°(eW^, e^U) and h\V, eU, eV, e^) = {W, eU, eW, e^U), 
we can write (64) as 

b^{eW, e^U)dtW =-c^diW + b\W, eU, eW, e'^U)diW 
e 

(99) + J^o{W, eU, eW, e^U) + eJ^i{W, eU, eW, e^U), 
where 

Fo{W,eU, eW, e^U) = /o(W^ + d^{W), eU, e{W + d^{W)),e^U) 

(100) - b^{eW, e^W)d^{W, eU, eW, e^U) + b^ {W, eU, eW)did<^iW), 
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and 

(101) J"i(W,eC/, eW, e^U)) = f^{W + di^{W),eU, e{W + d^{W)),e^U). 

In the following Proposition, the constant Csob is defined to be the e- 
independent in the weighted Sobolev inequality || • < Csobll • \\h^^ 

which holds for i > 3/2 + 1 and < e < eo (see Lemma A. 7 in [10] for a 
proof). 

Proposition 4.3. Suppose -1 < 6 < -1/2, eo > 0, s G No, > 0, 

Ki < Ko/i2^Csob), r > 2Ki/Csob, R > I6r+R, k > 3+s, a,w^ G H^^^, 

o o 

supp a C Br, 3^"^ G H^^^ , Ia^ G -f^l-i ■ Let ui , dtUe and wf be the initial 

o o o o 

data constructed in Proposition 4.1, which, by choosing < 1 small enough, 
satisfies 







for all e G (0,eo]. Then there exists a T > independent of e & (0,eo]; o^'^'d 
maps 

= Wg,6ui^,a,,wif G XT,,s,k,S-i < e < eo 

such that 

(i) T,>T forO<e< eo, 

(ii) We is the unique solution to (99) with initial data 

T 



\ o o o ° " 

(iii) 

— 1, £ -'-'(5 — 1,E 

and 

max{||eu^^'(t)||L^, \\ea,(t)\\L^ ,\\ew\t)U^} < 2Kq 

for all (t,e) G [0,T] x (0,eo], 
(iv) for each e G (0, eo], if 

limsup ||We(i)||^i,oo < oo , 

and 

sup {||eu*^'(t)||L^,||ea,(t)||L^,||ei(;*(t)||L^} < 2i^o 

0<t<T^ 
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then the solution W^{t) can he uniquely extended for some time T* > 
T 

(v) for any time which is strictly less than the maximal existence time 
and for which 

sup {||eu:^'(t)||L-, ||ea,(t)||L-, ||ew;*(t)||L-} < 27^0 

0<i<Tj 

holds, the support of satisfies 

supp a,{t) C Bji^ Vt G [0,1;] 

where := IGsupg^^^j, (t)||Loo + R, 

(vi) supp ae{t) C Bj^ for all {t,e) £ [0,T] x (0,eo], 

(vii) dtUe = e~^ul^, and djUe = + 5\5\di(^{ae) where ui = ui + 

(viii) the triple {ui ,ae-,wl} determines, via the formulas (7), (9), (14), and 
(15), a solution to the full Einstein- Euler system (1) in the harmonic 
gauge (26) on the spacetime region = x [0,r], and 

(ix) the conclusions (vii)- (viii) continue to hold on any region of the form 
D^=R^ X [0,Te] provided supp a,{t) C % for all < t < f^. 

Proof, (i)-(vii): First we observe, that proof of statements (i)-(vii) follow 
from a slight modification of the proof of Proposition 3.4 in [11]. 



(viii)-(ix): Let satisfy the initial value problem 

(102) v''dkA = : ^.(0) =X3if/2(0), 

and define 

Me = = V f (^ViV"- + - 



e 

Next, we observe that (52) contracted with yields 

2 

(1 - 2eMe)v^djMe - ^v^B.oMe + e(l - 2d\fe){x4R " ^)v'v'r^jVk = 0. 

Multiplying this equation by ip^ then gives 
(103) 

2 

(1 - 2eMe)vh>djK - ^^'9jaK + e(l - 2eMMxm - l)vWijVk = 0. 

From statements (i)-(vii), we have that 

(104) wiit) e XT,s,k,s-i 
and 

(105) \\w\t)\\H.<\\w\t)\\j,, <1 V(t,e)G[0,T)x(0,e]. 
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Since 

v^d^ = (1 + ewt)dt + widi, 

we get from (104), (105), and the hyperbolic equation (102) that ^ G 
XT,s,k,S-i and 

11^.(^)11//^ <i V(t,e)e[o,r)x(o,6]. 

From the finite propagation speed property of hyperbohc equations, we con- 
clude that there exists a time T* G (0, T) such that 

(106) Mt)\Bj, = 1, and suppVe(i) C B^j^ V {t,e) G [0,r,) x (0,eo]- 
In particular, this implies that 

(107) M,it,x) =K{t,x) V (t,x,e) G [0,r,) X X (0,eo], 
and 

(108) {x4Ri^)-l)A{t,x) = V(t,x,e) G [0,r,) xR3 X (0,eo]. 
Using (108), equation (103) reduces to 

2 

(1 - 2eJ\fe)vb>diK - ^v^d-juMe = 
2nh 

for all (t,x,e) G [0,r*) x x (0,eo]. But A4(0) = from the choice of 
initial data which implies that A''e(0) = 0. By the uniqueness of solutions to 
hyperbolic equations, we conclude that J\fe{t., x) = for all (t, x, e) G [0, T*) x 
R3x(0,eo], and hence TV, (t,x) = for ah {t,x,e) G [0, T*) x x (0, eo]. This 
implies that the fluid velocity normalization v'^Vi = — 1/e^ is satisfied for all 
(t, X, e) G [0, T^,) X X (0, eo]- Using this and the fact that supp ae{t) C i?R 
for all (t, e) G [0,T^,) x (0,eo], it is not difficult to verify from the evolution 
equation (51)-(52) that 

{v'^it^x) = l+ew^^{t,x),v\t,x) =wl{t,x),p,{t,x) := (4i^n(n+l))""a^"(t, x)} 

satisfy the Euler equations (46) (or equivalently (43)) for all {t,x,e) G 
[0,T^,) X X (0,eo]. With the Euler equations satisfied, the remainder 
of the proof follows as in the proof of Proposition 3.4 in [11]. □ 

5 Limit equations 

In this section, we describe the limit equations that govern the gravitational 
and matter fields in the limit e \ 0. We show in the next section that 
solutions to these equations approximate the solutions to the full Einstein- 
Euler equations up to a remainder term that is of order e. 
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5.1 Fluid limit equations 

The fluid limit equations are 

(109) dta = —w dja dju) , 

2n 

(110) dtw^ = --^d'^a - w^diw^ - XarQ^^, 

(111) Al- = p (p := {AKn{n + 1))-"q2") . 
Proposition 5.1. Let k, s, R, 5, a, and w be as in Proposition 4.3. Then 

o o 

there exists a maximal time Tq^ > and a unique solution 

a,w' e cO([o,ro^O,i^ti)nci([o,To),if|r/), 
!> G cO([o,To^^),i/^2^nC7H[o,ro^),/f^^), dt~^ e c\[^,t^),hI+1) 

to (109)-(lll) satisfying d(0) = a and w\^) = . Moreover, 

o o 

and 

supp d(t) C -B^(t) 
where R{t) = R + t supQ<j,<( ||t()^(s)||Loo . 

Proof. The proof follows from a trivial modification of the proof in Propo- 
sition 3.7 in [11]. □ 

Remark 5.2. Since R > R, it is clear from Proposition 5.1 and the weighted 
Sobolev inequality (see Lemma A. 7 in [10]) that there exists a time Tq G 
(0, Tm) such that 

(112) supp d(t) C for aU t G [0,ro]. 
In particular, this shows that 

(113) p{t)X4R = pit) for all t G [0, To], 

which in turn implies that the pair {p{t),w^ (t)} satisfies the Poisson-Euler 
equations (2)-(4) on the time interval [0,To]. 

5.2 Gravitational limit equations 

The gravitational limit equations are defined by 

(114) 9fX, = -^C'diX, + (0, 0, xlfj^, 
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where 



e=0 



j^, = ixl,xl,xi^] 

Proposition 5.3. Let 5, k, a, dtuf , djUe be as in Proposition 4.3, and 

o o o 

(115) X,(0) = (edtui^, diui^ - diA-\6i6ip),0] 

\ o o o / 

Then there exists a unique solution 

X, G C°([0, oo), ifti) n C\[0, oo), Htl) 

to (114) with initial data (115) that satisfies G -'^oo.s.fc.^-i ^^^^^ the esti- 
mates: 

(i) 



|X.(t)||H| , +e||5iX,(t)||^.-i <e^*, 
|X:^-(t)||^o.^ + ||Z)X:^-(t)||^.-,^ < e^*6, 



/or a// (t, e) G [0, oo) x (0, eo] and some fixed constant C > 0, and 
(ii) for any A > 0, 

(e + 1)3/2 

for all (t, e) G [0, oo) x (0, eo]. 



evil ~r i 



Proof. Since —1 < 5 < —1/2, it follows from Lemma A. 11 of [10] and 
Proposition 4.1 that 

(116) l|X.(0)||„. ^,,<||X,(0)||„. <1. 

o — 1 ' o — 1 

This inequality together with the weighted energy estimates (see Lemma 7. 1 
in [10]) gives 

(117) ||X,(t)||^. <e^*||X,(0)||.. <e^\ 

— 1 — 1, e 

for some fixed positive constant C. 

From the evolution equation (114) and the choice of initial data, we see 
that 

dt{diX'l{t)-djXfJ=0 and djX%{0)-djXf^^ = 0, 
which implies that 

(118) djXl{t)-djXlit) = 0. 
Also, it is not difficult to show that 

(119) dt{djXl^ -eX\{) = 
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follows from (114). Combining (118) and (119) then yields 

(120) diXl^{t) = e{xl{t)-Xl{0)). 
Next, we note that 

(121) WXl^h- < \\DXl^\U ^ +e\\Xl^h2 ^ 

follows from the weighted Sobolev inequalities (see Lemma A. 7 in [10]). 
Collecting the estimates (117), (120), and (121), we arrive at 

\\Xnt)\\Lr, + \\DXl^U <e^'e. 

"i^ — 1,6 

To prove the last two estimates for X^^^ and Xj^ , we observe that X^^^ 
and satisfy the wave equations 

(122) e^S^xjf^ - AXif^ = and e^^^x^f^ - AX)f^ = 0. 

Since the initial data for these equations satisfy (116) and —1 < 5 < —1/2, 
we can apply the weighted dispersive estimates from Theorem 1.1 in [4] to 
obtain 

, , .. , uUmu^, ^ + \\dixtmHri 

(123) \Dixyit,x)\ < ■ — , 

' (l + t/e + |x|)V|l + |i/e-k|| 



and 



(124) \Dixl{t,x)\< 



4,.r 



(1 + t/e + |x|)v^|l + |t/e- |x|| 
for < ^ < /c - 3/2. But 

1 + \t/e\ < 1 + \t/€- \x\ + \x\\ < 1 + A + |t/e - < (A + 1)(1 + \t/e-\x\ 
for \x\ < A, and so the inequalities (123) and (124) imply that 



for < ^ < /c - 3/2. □ 

Remark 5.4. From the initial value problem (114)-(115) (see (120)) and 

Propositions 4.1 and 5.1, it is not difficult to verify that 

(125) 

:=lx:i + (5)5^^3^^-2A-ia,3i-'4''5S^+'J4<54A"'(P+^lj/-')) + (^'(i)-^'(0))'^l<5i 

£ o 

satisfies the identities 

(126) = -^X'l + 6i6idt^, S/fi^J = Xy^^ + 5l6idi^, 
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and the wave equation 

(127) e^df^i^ - Au:^' = -5l6ip + dldie^d^^, 
with initial conditions 

(128) a^^l^^o = ^i^W - 2A~%ii'5^:6^] + 6i6im^^^+A~'dW), 
(129) 



5.3 The combined system 

Collecting the fluid and gravitation limit variables into a single vector 

(130) Ye = (Xe,a,^^^^S^)^, 

we can write the combined gravitational-fluid limit equations in the following 
form: 

(131) dtY, = ^c^diY + bidiY, + Jo + c^diu 
where 

(133) 

u = (cu^^o;^, 0,0,0,0)^, u;'^ = dM6l u:\' = diA-'{2pw'6^]5i^), 
and 

(134) -Fo =( - 2pw'^6^i5i\ -dim6l6^,xi^ ,0, -x^Rd'^, of. 



6 The fast Newtonian limit 

We begin by deflning the error between the limit and the full solution 

by 

(135) We = Y, + e{uj, + Z,). 
Next, we let 

(136) T = ToiYe,0,0,0), and 6^ = 6^(y,, 0, 0, 0), 
and observe that 

(137) T = To + eTi, and V = bi + eb{, 
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where 



(138) bi='-r 0), 4Xf^ 

^ ^ ^ \ Oy 

and 

(139) J-i = - (4(xf^9|j4> - /57?,,X:^) JiJi, 0, 0, 0, P.f. 



with 
and 



•^1 = -X4fl(^4,t + ^KlX^J') 



Ml 
4 > 



Now, let T* = min{To,T} where T and Tq are as defined in Propositions 
4.3 and Remark 5.2, respectively. Then by Propositions 4.3, 5.1, and 5.2, 
for any e € (0, eo], the error Zf:{t) (see (99), (131), and (135)-(137)) satisfies 
the initial value problem 

(140) b%z, = ^c^diz, + bidiz, + 

(141) Z,(0) = -{W,{0) - YM) - ^(0), 

e 

on the interval < t < T^,, where 

(142) 6° = 60(et^e, e^U,), hi = b'{W,, eU,, eW,, e^U,), 

(143) = J-o(T^e, eC/e, 6^^., e'C^e) + eJ'iiW,, eU,, eW,, e^U,), 

and 
(144) 



= bidiio - b%Lo + ^-^ — 9/y + — — + ' ^ edtY + b{diY + Tx . 



e e 

Proposition 6.1. Let 6, k, s > 2, T, and VFe(t) be as in Proposition 4.3, 
To as in Remark 5.2, Y^{t) as defined by (130), and = min{To,T}. Then 
for eo > small enough 

||VF,(t)-y,(t)|| .-2 <e 

o — 1 ,e 

for all (t,e) G [0,r=,) x (0,eo]. 

Proof. By Propositions 4.1, 4.3, and 5.1, there exists a positive constant Cq 
such that 

(145) \\ZM\\h^ <Co forahee (0,eo]. 
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Next, choosing eg small enough, it follows directly from Propositions 3.5 and 
3.6 of [11] and Propositions 5.1 and 5.3 of the previous section that 

(146) bidiu - h%u + ^l^diY, + - ^0 



< 1 + 11^. 



,5-1, E 

fc-2 



e e 

for all (t,e) G [0,r*) x (0, eo] provided \\Z,{t)\\ k-2 < 2Co/e. Also, from 

(5 — 1 , e 

Lemmas A.l and A. 4 of [10], and Proposition 5.1 and 5.3, we see that 



(147) 

for all (t,e) G [0,r,) x (0,eo]. 
Defining the energy norm 

III • ||U-2,<5,e '■ = 

we see via Proposition 4.3 that 

< I 



<i + 



\a\<k-2 



\\k-2,5-l,e ^ 



uniformly for (t, e) G [0, T*) x (0, eo]- Setting = eZ^, the evolution equa- 
tion (140) and the weighted energy estimates (see the proof of Theorem B.l 
in [11]) in conjunction with Proposition 4.3 and the estimates (146)-(147) 
show that there exists a fixed constant Ci > such that 

-lZ,{t)\\k-2,5-l,e < Cl |||^.(t)||U-2,5-l,. + e+ ( ^ 



for all t such that \\Zf 
then show that 



\\Z. 



k~2,5~i,e < e'-^'Coe + e / e 





k~2,5~i,e ^ 2Co. Gronwall's inequality and (145) 

ds 



* pCi{t-s) 



< ee 



Co + 1 + 



1 



(e + s)3/2 



(e + s)3/2 
ds 



<ee^i*(C7o + 3), 



again for all t such that ||-2e(t)||fc-2,5-i,e < 2Co. Therefore choosing eo > 
small enough we obtain 



for all (t, e) G [0, Xi<) x (0, eg], and the proof is complete. 



□ 



We are now ready to prove the main theorem. 
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Proof of Theorem 1.1. 

(i) : Since the ADM mass is conserved, statement (i) follows directly from 
Corollary 4.2 and Proposition 4.3. 

(ii) -(iv): From the definition of and 1^, we have 

\\a,it)-ait)\\j^,^2 +\\wiit)-w'it)\\^ < \\W,{t)-Y,{t)\\j^,-2^, 

and hence, by Proposition 6.1, 

(148) \\a,{t) - ait)\\^,-2 + \\wi{t) - mn'^^^ + llw^^WII^^-^ <e 

— 1.6 — 1, e — 1, e 

for all (t,e) G [0, T*) x (0, eo]. Also by the weighted multiplication Lemma 
(see Lemma A. 8 in [10]) and Propositions 4.3 and 5.1, we have 

(149) \\p,it) - /5,(i)|| „.-2 < ||ae(t) - ait)\\Ht' 

o — 1,€ o — l,e 

for all te[0,T^)x (0,eo], while 

wpeit) - m\\H>^-^ + - w\t)\\Hk-2 + \\wfmH>^-^ 

(150) < \\p,it) - m\\Ht' + - + w^tmHt' 

— 1,6 — 1, e — 1,6 

is a consequence of Lemma A. 11 and equation (A. 24) of [10]. Combining 
the inequalities (148)-(150), we arrive at 

(151) \\p,{t) - mWm-^ + Hit) - w'it)\\jj.-2 + H{t)\\H,-2 < e 
for all te[0,n)x (0,eo]. 

Next, we observe that 

\\di^S)-dim\\^,-2 +||e5tcD,(t)|| <e 

(5 — 1,6 (5 — 1,6 

for all {t, e) £ [0, T,,) x (0, eo] by Propositions 4.3 and 6.1, and Lemmas 3.2 and 
3.3 of [11]. From the above estimate, the identities Wj-'^ = djUe + ^l^idi^e 
and U4 ^ = edtUe (see Proposition 4.3), and the relations (125)-(126), we get 

<\\wl{t) - xl{t)\U-2 +\\iCiJt)-Xl{t)\\^.-2 +e 

<\\W,{t)-YM\Ht' 
and hence, by Proposition 6.1 and Lemma A. 7 of [10], 

\\K'it)-K'mLl +\\djul^it)-dj^l^it)\\^.-2 +\\edtU^^^{t)-ediu^,^it)\\^.-2 <e 
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for all (t, e) G [0,T*) x (0, Eq]- Finally, it follows from Lemma A. 11 and 
equation (A. 24) of [10], and the above estimate that 

for all (t, e) G [0,T^) x (0, eq]. This completes the proof. □ 



References 

[1] R. Bartnik, Phase Space for the Einstein Equations, Comm. Anal. 
Geom. 13 (2005), 845-885. 

[2] U. Brauer and L. Karp, Local existence of classical solutions for the 
Einstein- Euler system using weighted Sobolev spaces of fractional order, 
C. R. Acad. Sci. Paris, Ser. I 345 (2007), 49-54. 

[3] G. Browning and H.O. Kreiss, Problems with different time scales for 
nonlinear partial differential equations, SIAM J. Appl. Math. 42 (1982), 
704-718. 

[4] P. D'Ancona, V. Georgiev, and H. Kubo, Weighted decay estimates for 
the wave equation, J. Differential Equations 177 (2001), 146-208. 

[5] H. Isozaki, Wave operators and the incompressible limit of the com- 
pressible Euler equations, Comm. Math. Phys. 110 (1987), 519-524. 

[6] S. Klainerman and A. Majda, Compressible and incompressible fluids. 
Comm. Pure Appl. Math. 35 (1982), 629-651. 

[7] H.O. Kreiss, Problems with different time scales for partial differential 
equations. Comm. Pure Appl. Math. 33 (1980), 399-439. 

[8] M. Lottermoser, A convergent post-Newtonian approximation for the 
constraints in general relativity, Ann. Inst. Henri Poincare 57 (1992), 
279-317. 

[9] T. Makino, "On a local existence theorem for the evolution equa- 
tion of gaseous stars", in Patterns and Waves, edited by T. Nishida, 
M. Mimura, and H. Fujii, North-Holland, Amsterdam, 1986. 

[10] T.A. Oliynyk, The Newtonian limit for perfect fluids. Comm. Math. 
Phys. 276 (2007), 131-188. 

[11] T.A. Oliynyk, Post-Newtonian expansions for perfect fluids. Comm. 
Math. Phys. (accepted). 

[12] A.D. Rendall, The initial value problem for a class of general relativistic 
fluid bodies, J. Math. Phys. 33 (1992), 1047-1053. 

[13] A.D. Rendall, The Newtonian limit for asymptotically flat solutions of 
the Vlasov- Einstein system. Comm. Math. Phys. 163 (1994), 89-112. 

[14] S. Schochet, Symmetric hyperbolic systems with a large parameter. 
Comm. partial differential equations, 11 (1986), 1627-1651. 



T.A OLIYNYK 



27 



[15] S. Schochet, Asymptotics for symmetric hyperbolic systems with a large 
parameter, J. differential equations 75 (1988), 1-27. 

[16] S. Ukai, The incompressible limit and the initial layer of the compress- 
ible Euler equations, J. Math. Kyoto Univ. 26 (1986), 323-331. 



